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We study how the multiplicity of the Fermi surface affects the zero-bias peak in conductance 
spectra of tunneling spectroscopy. As case studies, we consider models for organic superconductors 
re-(BEDT-TTF) 2 Cu(NCS) 2 and (TMTSF) 2 C10 4 . We find that multiplicity of the Fermi surfaces 
can lead to a splitting of the zero-bias conductance peak (ZBCP). We propose that the pres- 
ence/absence of the ZBCP splitting is used as a probe to distinguish the pairing symmetry in 
k-(BEDT-TTF) 2 Cu(NCS) 2 . 



I. INTRODUCTION 

An unambiguous determination of pairing symmetry in 
unconventional superconductors is crucial to understand 
the pairing mechanism of superconductivity Strong ev- 
idences suggesting d x 2_ y 2-w&ve pairing symmetry in the 
high-T c cuprates have been provided using several phase- 
sensitive probe a 1 ! 2 ' 3 including tunneling spectroscopy via 
Andreev surface bound states (ABS's)4*ii*§*£ The tun- 
neling spectroscopy via ABS's enables us to detect the 
sign change in the pair potential as well as its nodal 
structureA^ This state, which originates from the in- 
terference effect in the effective pair potential of the 
eLj,2_ a 2-wave symmetry through reflection at a surface or 
an interface, have significant influences on several charge 
transport propertiesi2iiiii2. The existence of ABS's, 
which manifests itself as a distinct conductance peak 
at zero-bias in the tunneling spectrum (zero-bias con- 
ductance peak, referred to as ZBCP), has been actually 
observed not only in the high-T c cuprates^ but also in 
ruthenates 1 ^* 1 ^, heavy fermion systems 1 ^, and more re- 
cently MgCNig, 1 - In this context, it is of great interest to 
investigate whether the ZBCP due to the ABS's can be 
observed in organic superconductors 1 ^ such as k-(BEDT- 
TTF) 2 AT and (TMTSF) 2 A^a. 

The tunneling spectroscopy via ABS's can be used to 
determine the pairing symmetry if one can prepare well- 
treated surfaces with arbitrary orientations in the su- 
perconducting plane. For high-T c cuprates, which has 
a d x i_ y i pair potential, it is theoretically shown that the 
ZBCP should be observed most prominently for (110) 
surfaces or interfaces. Moreover, it has been clarified that 
the ZBCP may be observed due to atomic-scale rough- 
ness even in the (100) surfaces i2L22*22i In fact, Iguchi 
et al. 24 » have measured the ZBCP for Ag/YBCO ramp- 
edge junctions with various orientations, where the in- 
jection direction varies continuously from (100) to (110) 
interfaces. The height of the ZBCP has shown to vary 
according as the misorientation angle from a-axis within 
the plane. 

As regards organic superconductors such as k-(BEDT- 



TTF^X, the pairing symmetry of the pair potentials 
still remains to be a controversial problem. It has in- 
deed become an issue of great interest whether k-(BEDT- 
TTF^X has a d-wave pair potential similar to high- 
Tc cuprates. There is now a body of accumulating ex- 
perimental evidences suggesting that k-(BEDT-TTF)2A" 
have anisotropy in the pair potentialiS^SSiSL 2 ^ 2 ^ 3 ^ Ear- 
lier theories support d x 2_ y 2-w&ve pairing ^i* 3 ^ 3 . while a 
recent thermal conductivity measurement suggests d xy - 
wave pairing^ Concerning this issue, two of the present 
authors have theoretically shown that a d^y-like pair- 
ing may slightly dominate over d x 2_ y 2 pairing when 
the dimerization of the BEDT-TTF molecules is not so 
strong. 35 According to previous studies^Sii if the pair- 
ing symmetry of k-(BEDT-TTF) 2 X is d-wave, ABS is 
expected to be created at surfaces for arbitrary injec- 
tion orientations. However, a scanning tunneling mi- 
croscopy (STM) experiment for k-(BEDT-TTF) 2 A" by 
Arai et aliSi showed the absence of ZBCP for arbitrary 
injection angle from the c-axis in the &c-plane, which is 
in contrast with the case of the high-T c cuprates. The 
presence/absence of the ZBCP of d-wave superconduc- 
tors is sensitive to several factors: (i) roughness effect 
of surfaces or interfaces, (ii) random impurity scatter- 
ing effect near the interfaces, (iii) the shape of the Fermi 
surface, and (iv) the degradation of surfaces. The dis- 
appearance of the ZBCP in d-wave superconductors due 
to reason (i) has been studied previously2L2£. Depend- 
ing on the shape of the Fermi surface and the geome- 
try of the surface, the atomic-size wave nature of the 
zero energy ABS's (ZES), i.e., the oscillatory behavior 
of the wave function of ZES induces an interference ef- 
fect which locally destroys the ZBCP. In fact, it is by 
no means easy to make well-oriented cleavage surfaces in 
organic materials, so this point may be important. As 
regards point (ii), Asano et al^ have shown, both from 
analytical and numerical calculation beyond quasiclassi- 
cal approximations 3 ^, that impurity scattering near the 
interface in the high-Tc cuprates can induce a splitting 
or a disappearance of the ZBCP. 

As for point (iii), we have recently studied the disap- 
pearance of ZBCP due to the warping of quasi- ID Fermi 
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surface as in (TMTSF) 2 A.22i The results indicate that 
the ABS's are sensitive to the shape of the Fermi sur- 
face. However, most of the theoretical studies on tun- 
neling spectroscopy via ABS's up to date have been per- 
formed for single band systems. It has not been clarified 
how the multiplicity of the Fermi surface influences the 
ZBCP. 

Motivated by this point, here we investigate the 
surface density of states in systems having multi- 
ple Fermi surfaces, where we focus on two organic 
superconductors as case studies, namely, k-(BEDT- 
TTF) 2 Cu(NCS) 2 and (TMTSF) 2 C10 4 . The Fermi sur- 
face of k-(BEDT-TTF) 2 Cu(NCS) 2 , which has been de- 
termined by Shubnikov-de Haas experiment^ consists 
of two portions separated by small gaps. The Fermi sur- 
face of (TMTSF) 2 C1C>4 is also separated by a small gap, 
which is due to anion ordering. In this paper, we extend 
our previous studies^ on anisotropic triangular lattice 
by taking into account these multiplicity of the Fermi 
surface. 

The organization of the paper is as follows. The 
formulation of calculating the tunneling spectrum on 
anisotropic triangular lattice is presented in Sec. [n] In 
Sec. IIIII results of the numerical calculations are dis- 
cussed in detail. Finally, we summarize the paper in 
Sec. HV| 



II. FORMULATION 




(b) d x 2. y 2-wave d X y -like-wave 




FIG. 1: (a) Schematic of (100) surface in xy-p\&ne with next- 
nearest neighbor hopping t' . (b) Cooper pairs in real space 
for d x 2_ y 2 and d xy - like- wave pairings. 



In the present study, we start from an extended Hub- 



bard model given by 

+ ^2( £ i ~ ^) c \,a c i,<yi (!) 

i.a 

where ct creates a hole with spin a —\,\ at site 
i = {i x ,i y ). As a model for k-(BEDT-TTF) 2 Cu(NCS) 2 , 
each site corresponds to BEDT-TTF molecule dimmers. 
We consider five kinds of hopping integrals, t x (= t), 
t x ' , ty, ty> j and t' in the xy plane on the anisotropic 
triangular lattice as shown in Fig. da). In order to 
reproduce the shape of Fermi surface for k-(BEDT- 
TTF) 2 Cu(NCS) 2 and (TMTSF) 2 C10 4 , we adopt the val- 
ues of (i) t' = Q.8t, t y ' = t x , t x > = ty^L and (ii) 
t y = O.lt, t x ' — t x , tyi = £ y > 41 i 42 i 43 respectively. Two 
subchains in the x direction alternatively have the site 
energy Si = +E g , —E g , +E g , ... in the y direction. 44 ! 45 
The chemical potential /i is determined such that the 
band in k-(BEDT-TTF) 2 Cu(NCS) 2 [(TMTSF) 2 C10 4 ] is 
half- filled [quarter- filled] . The effective attraction V is 
assumed to act on a pair of electrons. 

By solving the mean-field equation for a unit cell with 
A ? l(= 500) sites in the x direction and two sites in the 
y direction, we obtain the eigenenergy E v . In terms of 
the eigenenergy E v and the wave functions u\, v%, the 
Bogoliubov-de Gennes equation for the (100) surface in 
the xy plane is given by 

with 

Hij(ky) = - t x rj + S jxtix+ i - t x >rj^5 jxjix+ i 

+ {(-iyy+i Eg -n}5 ix jJ iv , jy , (3) 
where we define rj + = |{1 + (— 1)^+ 1 h} and r\_ = i{l — 

As for plausible pairing symmetries in k-(BEDT- 
TTF) 2 Cu(NCS) 2 , we consider d x 2_ y 2-wave pairing given 



3 



(a) d x 2.y2-wave 
ii — ■ — ■ — ■ — i— 



d-wave 







t X '/t x =2 






y/t,=as 









(b) d X y -like-wave 




o 

k c /lt 





k,.l% 



FIG. 2: The Fermi surface and d-wave pairings, (a) d x 2_ y 2 
and (b) d xy - like- wave. 



by, 



FijQty) = A x r] + 5j x , ix+1 + Aaj/jy-^.^+i 



A y e 



A^e 



A y e" 



(4) 



and d^-like pairing given by 



Fij(k y ) = A x r] + Sj xjix+ i + A x ,rj-6j xt i x+ i 
+ A y e- 2ik * 5 ^ V+ 6 jv , iv+1 

-aA p e- 2ik ^8 jx , ix+1 S jyiiy+1 
+ A x r)-S jxiix -i + Axirj+Sj^-! 
+ A v e* k ^r 1 _6 jv , iv _ l 
+ Ay^^n+K^-x 



aA p e 2ik ^-'6 jx , ix 



-l°jy,iy- 



(5) 



with a — 0.8f in accord with Ref. 135L The pairing in real 
space is shown in Fig. db). Here, we select A x i = A y 
and A x — A y > = A p = Aq, where Ao is a bulk value. 

The upper and lower panels of Fig. [21 show the d x 2_ y 2 
and ci^-like pair potentials in momentum space along 
with the Fermi surfaces. For t x > = t x , the Fermi sur- 
face is continuously connected. In the actual k-(BEDT- 
TTF) 2 Cu(NCS) 2 , however, BEDT-TTF dimmers are 
further dimerized so that t x i =fi t x , which leads to 




(b) t'=-0.08t 
1 




FIG. 3: (a) Cooper pairs with d-wave symmetry separated by 
two lattice space and (b) quasi-lD Fermi surface in t y — 0.1t x 
and t' = -0.08^. 



a splitting of the Fermi surface at around (fc c , kf,) — 

(±7T,±7r/2). 

In (TMTSF) 2 C10 4 , the orientational order of the an- 
ions CIO4 doubles the unit cell, leading again to a split- 
ting of the Fermi surfaced Although the pairing symme- 
try for (TMTSF) 2 X remains to be undetermined, here 
we assume singlet d-wave pairing as an example in which 
the multiband effect is prominent. In this case <i-wave is 
a pairing separated by two lattice spacings given by 



F, 



A 2x S 



j x i ^x ~1~2 



A 2x 8 ]x 



(6) 



The anion potential E g is estimated from experimental 
measurement of angle dependent magnetoresistance4£ 

In order to compare our theory with STM experiments, 
we assume that the STM tip is metallic with a constant 
density of states, and that tunneling occurs only to the 
site nearest to the tip. This has been shown to be valid 
through the study of tunneling conductance of uncon- 
ventional superconductors^. The tunneling conductance 
spectrum is then given at low temperatures by the nor- 
malized surface density of states^, 



P(E) = 



dwps(^)sech 2 



U3 + E 
2k B T 



dwpN(^)sech 2 



w - 2A 



2k B T 

psH - E [K l 2 ^ - E ») + Kl 2 % + E »)] 

kb,v 



(7) 
(8) 
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FIG. 4: Tunneling spectrum for (a) d x 2_ y 2 and (b) d xy -like- 
waves fixed in A x i — A x . 



Here ps(w) denotes the surface density of states for the 
superconducting state while p^(lo) the bulk density of 
states in the normal state. 



III. RESULTS OF IN-PLANE TUNNELING 
SPECTRUM 

In this section, we present the calculation results. 
First, let us focus on the model for k-(BEDT-TTF) 2 A. 
We examine the case of the tunneling spectrum at (100) 
surface on the xy plane as shown in Fig. ^ As seen in 
Fig- E] in the case of t x > = t x , where the Fermi surface is 
elliptical but continuous, there exists a distinct peak at 
zero energy, which resembles those obtained in previous 
theories assuming round shape Fermi surface. The ZEP 
arises because incident and reflected (including oblique 
incidence) quasiparticles normal to the surface feel op- 
posite signs of the pair potential, which results in a for- 
mation of the ABS. If we turn on the multiband effect 
by letting t x > ^ t x , the ZEP is found to split into two. 




FIG. 5: The t x > /t x v.s. the zero-energy peak splitting width 
5. 



This is reminiscent of the ZEP splitting originating from 
broken time reversal symmetry states A^^LiiSi 
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FIG. 6: Tunneling spectrum for variour A x i/A x in t x , jt x 
0.9. 



We have further studied the t x i / t x dependence of the 
ZEP splitting. In Fig. the width of the ZEP splitting 
S is plotted as functions of t x >/t x for d x 2_ y 2 and d xy -\ike 
pairings. S for the d xy -\\ks pair potential is almost pro- 
portional to t x //t x , and larger than that for d x 2_ y 2. In 
the regime of t x * jt x > 0.9, in particular, we see no split- 
ting for the d x 2_ y 2 pairing. Since t' x /t x is estimated to be 
~ 0.9jiS we may be able to distinguish between d x 2_ y 2 
and eicy-like pairings through the presence/absence of 
ZEP splitting. 

We have also performed similar calculation by letting 
A x > / A x deviate from unity, which should be the case 
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when t x deviates from t x . The results are plotted in 
Fig- EH for various A x i/A x with t x >/t x fixed at 0.9. In 
this case, we observe an overall shift of the splitted ZEP, 
while the magnitude of the splitting remains unchanged. 

Let us now move on to the model for (TMTSF) 2 C10 4 . 
In this model, as the Fermi surface becomes asymmetric 
with respect to k x — ► — k x transformation, some injected 
and reflected quasiparticles feel different signs and the 
ZEP appears in the tunneling spectrum^. When E g is 
turned on, a minigap opens at k y — ±7r/4. This effect 
again leads to ZEP splitting, of which the magnitude 
increases as E g is increased. [Fig. [7] 

Although it is by no means easy to pinpoint the origin 
of the ZEP splitting analytically, it can be qualitatively 
explained as follows. For single band models, the ZEP 
appears due to the sign change of the pair potential felt 
by quasiparticles at the interface. In multiband systems, 
injected and reflected quasiparticles have different band 
indices, so that an additional phase factor is expected 
to appear due to interband scattering. This additional 
phase factor induces the ZEP splitting as in the case of 
pair potentials with broken time reversal symmetry. 
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FIG. 7: Tunneling spectrum in t y = 0.1t x and t' = —0.08^ 
with anion potential E g . 



IV. CONCLUSIONS 



To summarize, we have investigated the multiband ef- 
fect on tunneling spectroscopy. As case studies, we have 
focused on models for k-(BEDT-TTF) 2 Cu(NCS) 2 and 
(TMTSF) 2 C10 4 . We find that the multiplicity of the 
Fermi surface can lead to a splitting of the ZEP. As 
regards k-(BEDT-TTF) 2 Cu(NCS) 2 , since t'Jt x is esti- 
mated to be ~ 0.9^ we can distinguish between d x 2_ y 2 
and c^y-like pairings through the presence/absence of 
ZEP splitting. 

As mentioned in the Introduction, however, a scanning 
tunneling measurement*^ actually find no ZBCP in the 
tunneling spectrum of k-(BEDT-TTF) 2 Cu(NCS) 2 . Since 
many experiments suggest the existence of nodes in the 
pair potential in k-(BEDT-TTF) 2 A, we believe that the 
absence of ZBCP is not because the pairing symmetry is 
a simple s-wave, but because of the roughness of the sur- 
face or the random scattering effect by impurities near 
the interface^, namely, point (i) or (ii) mentioned in the 
Introduction. As regards the roughness of the surface, 
we believe it is necessary to study tunneling spectroscopy 
of k-(BEDT-TTF) 2 Cu(NCS) 2 in the presence of atomic 
scale roughness as done by Tanuma et al. on a lattice 
model^i*2£ As for the issue of random scattering effect 
by impurities, Asano et al& have shown, both from an- 
alytical and numerical calculation beyond quasiclassical 
approximations, that impurity scattering near the inter- 
face in the high-Tc cuprates can induce a splitting or 
a disappearance of the ZBCP. From this viewpoint, it 
would also be interesting to study the impurity scatter- 
ing effect in k-(BEDT-TTF) 2 Cu(NCS) 2 . 

In order to clearly determine the pairing sym- 
metry, other complementary probes should also be 
used. Recently, we have shown that magnetotunnel- 
ing spectroscopy is a promising method to identify 
the detailed paring symmetry of the unconventional 
superconductors22i2L£i. It would also be interesting to 
apply this probe to k-(BEDT-TTF) 2 Cu(NCS) 2 . 

It is well known that ABS's have serious influence on 
the Josephson current. There are many works on Joseph- 
son effect in unconventional superconductors both from 
theoretical and experimental view point 52 . It is also a 
future problem to study Josephson effect in k-(BEDT- 
TTF) 2 Cu(NCS) 2 and (TMTSF) 2 C10 4 . 
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